The Optical Excitation of Zigzag Carbon Nanotubes with Photons Guided in 

Nanofibers 
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We consider the excitation of electrons in semiconducting carbon nanotubes by photons from 
the evanescent field created by a subwavelength-diameter optical fiber. The strongly changing 
evanescent field of such nanofibers requires dropping the dipole approximation. We show that 
this leads to novel effects, especially a high dependence of the photon absorption on the relative 
orientation and geometry of the nanotube-nanofiber setup in the optical and near infrared domain. 
In particular, we calculate photon absorption probabilities for a straight nanotube and nanofiber 
depending on their relative angle. Nanotubes orthogonal to the fiber are found to perform much 
better than parallel nanotubes when they are short. As the nanotube gets longer the absorption of 
parallel nanotubes is found to exceed the orthogonal nanotubes and approach 100% for extremely 
long nanotubes. In addition, we show that if the nanotube is wrapped around the fiber in an 
appropriate way the absorption is enhanced. We find that optical and near infrared photons could 
be converted to excitations with efficiencies that may exceed 90%. This may provide opportunities 
for future photodetectors and we discuss possible setups. 

PACS numbers: 78.67.Ch, 78.40.Ri, 73.22.-f, 78.30.Na 



I. INTRODUCTION 

The unique physical properties of carbon nanotubes 
and the flexibility they provide in selecting their char- 
acteristics offers great potential for nanotechnologj*i~— . 
Carbon nanotubes can be either semi-conducting or 
metallic, depending on their diameter and helical config- 
uration. They typically have nanometer sized diameters 
and a length of a few microns, although centimetre long 
nanotubes have been produced recently^. This makes 
them ideal ID systems that possess a ballistic conduct- 
ing channel^, no backward scattering, and energy levels 
that can be adjusted with external fields^—. Supercon- 
ductivity has also been observed in multi-walled carbon 
nanotubes and single carbon nanotubes have exhibited a 
superconducting proximity effect^—. Their applications 
range from extremely strong fibers and organic electron- 




ics^ to electrochemical sensor o 14 i 15 and photon detec- 



tor; 
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Carbon nanotubes are a form of carbon formed by 
rolling up a sheet of graphene into a cylindrical tube. 
An illustration of this is given in Fig. [TJ Here we fo- 
cus on the optical properties of carbon nanotubes. For a 
straight nanotube, inside a weak uniform classical plane 
wave field, these properties have been extensively stud- 
iedi£~— . Their quasi one-dimensionality means that their 
density of states exhibit Van Hove singularities and these 
contribute to strong optical absorption peaks. However, 
these results apply the dipole approximation, where it is 
assumed that the field does not vary along the nanotube's 
length. We extend this treatment by allowing for the spa- 
tial dependence of the field. This situation is particularly 
relevant when the electrons are delocalized in a tightly- 
confined field, such that the field varies greatly over a few 
hundred nanometres. The degree to which the electrons 



FIG. 1. (color online) a) Graphene lattice with the unit cells 
vectors labeled ai and a2. These vectors have the length 
a. Here the atoms in the A sublattice are red (dark grey) 
and the B sublattice is cyan (light grey). a c is the distance 
between neighboring atoms. The unwrapped unit cell for a 
zigzag (3,0) nanotube is shown shaded and the C vector de- 
fines the nanotubes circumference, b) A section from a zigzag 
(7,0) nanotube is shown with its unit cell shaded. T is the 
tangential unit vector of the nanotube's unit cell. 



are delocalized is a topic of ongoing research and var- 
ious studies have been done on the coherence length in 
nanotubes. Their results range from 10 nm to several mi- 
crons suggesting that the spatial field dependence is cer- 
tainly important for confined fields and may also be rele- 
vant for plane waves 2 ^— . The systems we are primarily 
interested in arc subwavelength-diameter optical fibers 
coupled to carbon nanotubes. The electrical field of a 
nanofiber is tightly confined and primarily exists outside 
of the fiber, in a large evanescent field^. Due to the pres- 
ence of a strong field in a relativity small volume, these 
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nanofibers arc ideal candidates to achieve a high optical 
absorption in atomic systems. For example, their interac- 
tion with atom-arrays has been studied 3 ^ - — . However, 
in contrast to such atom-fiber systems, the dipole ap- 
proximation can not be applied in the case of nanotubes 
since the optical field typically changes rapidly along a 
nanotube's length. In this paper we calculate the (inter- 
nal) quantum efficiency, i.e. the probability that a nan- 
otube, placed inside the evanescent field of a nanofibcr, 
absorbs a photon. Calculations for the external quantum 
efficiency, i.e. the efficiency for detecting the excitation 
with the photocurrent, are beyond the scope of this pa- 
per. However, it should be noted that important effects 
that could aid in this procedure, such as the avalanche 
effect, have been observed in carbon nanotubes 3 - 7 -. We fo- 
cus specifically on the example of zigzag nanotubes (see 
Fig. [IJ because they can be direct semiconductors. How- 
ever, our results are still representative of other semicon- 
ducting nanotubc types such as chiral nanotubes. We 
find that the absorption is extremely dependent on the 
nanotube's orientation. These results are highly relevant 
for the interface between any future nanoscale photon- 
ics and carbon nanotubes. If the absorption process is 
coherent the system may also be suitable as a quantum 
memory, that maps a photonic quantum state on to a 
coherent excitation of the nanotubc. 

Wc will be using the band-to-band tight binding tran- 
sition model for the carbon nanotube. This has proven 
itself to be very effective in determining the basic opti- 
cal properties of nanotubes but does not include effects 
due to excitons 3 ^ and electron-electron interactions^. 
Such effects give measurable corrections and there have 
been a few studies considering the exciton absorption 
strengt h 17 ! 23 ' 40 ! 41 . Nevertheless, the band-to-band model 
is suitable to determine the main contributions to optical 
absorption. 

This paper is organized as follows. We begin by giving 
an overview of nanotube properties and the calculation 
of their band structure in Sec. [TTJ Based on this we then 
evaluate the photon absorption by zigzag carbon nan- 
otubes in Sec. Mil In the setups considered in this pa- 
per the nanotubes experience fields that change strongly 
along their length, i.e. to calculate photon absorption 
we can not rely on the dipole approximation. We ob- 
tain general expressions for the absorption probability 
which arc then applied to cylindrical vacuum cladded sil- 
ica nanofibers in Sec. IIV1 and discuss different geometri- 
cal setups of nanofiber and nanotube. Possible photode- 
tectors that use these setups are then presented in Sec.lVl 
Finally, in Sec. IVI1 we summarize our results. 



II. THE TIGHT-BINDING MODEL FOR THE 
CARBON NANOTUBE 

Here wc will review the basic properties of carbon nan- 
otubes for completeness and layout the notation that 
we use in later sections. A single walled carbon nan- 



otube (SWCNT) can be thought of as a sheet of graphene 
wrapped into a tube, so we will start by describing the 
tight-binding model of graphene 4 ^. Graphene is a reg- 
ular 2D hexagonal Bravais lattice of carbon atoms and 
its structure is shown in Fig. [TJ We label the unit vec- 
tors of the graphene lattice ai and a.2. The length of 
these vectors is the lattice constant a which is related 
to the distance between neighboring carbon atoms, a c , 
by a = \/Sa c ~ 0.246 nm. Within each unit cell there 
are two carbon atoms, that we label to form the A and 
B sublattices. We can then define the unit vectors of 
the reciprocal lattice as bi and b2, with a^.b^ = 27r<5y. 
The first Brillouin zone given by these is also hexagonal. 
It has a selection of points with high-symmetry; one at 
the center, the midpoints of the hexagonal edges and two 
incquivalent types of corners. 

The well-established tight-binding model assumes that 
the electrons are tightly bound to the individual carbon 
atoms and the localized atomic orbitals are used as a ba- 
sis for expanding the wavefunction. Wc consider orbitals 
that contribute to states that lie within an optical range 
of energies around the Fermi level. These are the states 
that give the main contributions to the optical properties 
of the nanotube. Every carbon atom has four valence or- 
bitals (2s, 2p x , 2p y and 2p z ) that could lie in this energy 
range. For 2D graphene the (s, p x , p y ) orbitals combine 
to form hybridized sp 2 orbitals. These give the strong 
covalcnt bonds; primarily responsible for the binding en- 
ergy and clastic properties of the nanotubes. In the tight- 
binding model they result in <j and a* bands. However, 
their energy levels are far away from the Fermi level and 
hence they do not play a key role in the optical prop- 
erties that we are interested in. That role is played by 
delocalizcd ir and ir* bands that arc formed from the p z 
orbitals 1 -. Hence, we can ignore the a electrons and re- 
strict the tight-binding model to the 7r electrons. The 
Hamiltonian for this system is 

H a = -joY,(*Wj+ h - c -), (!) 

ij 

where —70 is the hopping amplitude and ij refers to near- 
est neighbors. Here, a\ and /3j are the creation operators 
for electrons in sublattice A and B, respectively. In this 
Hamiltonian wc have removed the constant energy contri- 
bution that corresponds to the Fermi level. We expand 
the wavefunction in terms of p z orbitals at every atom 
site and split this expression into two parts; one for each 
sublattice. The wavefunction for each state is then 

*(k,r) =]TC A (r A ,k) Pz (r-r A ) (2) 
+ ^C B (r B ,k)p z (r-r B ) (3) 

t'i! 

with rA,rB labeling the atom locations in sublattice A 
and B, respectively. The individual p z orbitals are given 
by the normalized wavefunctions p z (r) and each one has 
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a coefficient, represented with Ca and Cb- By using 
translational symmetry we can represent this as 



*(k,r) = c A (k)pf(k,v) + c B (k)pf (k,r) 
where the Bloch functions, pf and pf , are 
1 



#(k,r) = 



pf(k,r) = 



V N cells 

1 

V N ce u s 



^e ikr > z (r-r A ) 



(4) 

(5) 
(6) 



Here N ce iis is the number of unit cells in the graphcnc 
sheet. These bloch functions have the coefficients ca and 
cb- The states are labeled by their crystal momentum 
vector k. We now solve the time-independent single- 
particle Schrodinger equation 



#o*(k,r) = £(k)*(k,r). 
We define the quantity 



(7) 



(8) 



where q are vectors from an atom in the A sublattice to 
its neighboring atoms in the B lattice. This gives us the 
following quantities 



H AA = H BB = (P*\H \p*) = (p?|#o|p?) = 



H AB = H* BA = (pf\H \p^ 



S A 



-70</>k 



B 



S 



(pf\p B " 



BA 



M0 k . 



(9) 



Now the variational Schrodinger equation in matrix form 
is 

&£:)(2)-«w( s y")(s 

' (10) 

The above matrix equation can be solved to give the en- 
ergy of each state as 



E±(k) = 



±7o 



(11) 



where 



|0K 



1+4 cos 




4 cos 



(12) 



A typical value for 70 is 2.89eV such that the tight- 
binding model corresponds with experiments 3 - ' 43 ' 44 . We 
will keep u in the equations but for all plots and numer- 
ical calculations we assume that u = 0, i.e. there is no 
orbital overlap. In Fig. [5] this 2D dispersion relation is 
plotted as a contour. In Eq. (fTTj) the signs refer to the 



two relevant bands, the conduction band and the valence 
band. The coefficients are found to be 



(13) 
(14) 
(15) 




«|0k|)' 



(16) 



Now we know the relevant band structure of graphene 
and their wavef unctions, we need to obtain the en- 
ergy states of the nanotubes. To do this we use the 
zone-folding approximation. This assumes the nanotubc 
bands are the same as graphene but with limited k, due to 
the ID nature of a carbon nanotube. There are a variety 
of ways available to wrap the sheet up into a nanotubc, 
each of which result in very different properties. The 
nanotubes are characterized by a vector in the graphene 
plane that corresponds to the circumference of the nan- 
otube and is called the chiral vector C = niai + nift-i 
(0 < |ri2| < tt-i) (see Fig. [TJ. It gives the relative position 
of two graphene atoms that become 'identical' when the 
graphcnc is rolled into a nanotube. We will use these 
parameters in the standard form (n 1 ,n 2 ) to label each 
type of nanotube. This immediately defines some ba- 
sic geometric properties such as the tube's circumference 
and radius R t = a^J n\ + nin 2 + n^fl-n. We also define 
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FIG. 2. (color online) Graphene Band structure and subbands 
for (7,0) nanotube bands. The transitions between bands al- 
lowed with the electric field perpendicular to the nanotube 
can occur between neighboring subbands. 

the translational vector, perpendicular to C, that corre- 
sponds to the direction along the tube, T = t\a.\ + £28-2- 



Using the greatest common divisor (gcd), we define 
ti = (2n 2 + nx)/N R , t 2 = -(2m + n 2 )/iV fl and JV fl = 
gcd(2ni + n 2 , 2n 2 + Hi). The two vectors, C and T, de- 
fine the unit cell of the nanotube. Within each nanotube 
unit cell there are N G = 2(n\ + nin,2 + tv^/Nr graphenc 
unit cells and, hence, N G = 2N G carbon atoms. In a 
nanotube of length L there are Nl = L/ |T| nanotube 
unit cells. For the nanotube's reciprocal lattice we define 
Ki = (tib 2 - t 2 b 1 )/N G and K 2 = (n 2 bi - ni h 2 )/N G 
such that K^.T = K 2 .C = and K^.C = K 2 .T = 2tt. 
These give the allowed vectors in the SWCNT's Brillouin 
zone to be a set of N G ID 'cutting lines' with values 



>'ir 



K2 

K 2 



(17) 



with fx = -N G /2 + 1, . . . , 0, . . . , N G /2 and -7r/ |T| < 
fen < 71"/ |T|. It is the periodic boundary condition along 
the circumferential direction of the tube that causes the 
wave vector to become quantized and each discrete cut- 
ting line is labeled by the azimuthal quantum number 
fx. For short nanotubes the wave vectors are also quan- 
tized along the nanotube's length causing discrete en- 
ergy levels to be formed^. These discrete values have 
fc|| = 2nj/L — it/ |T|, for an integer j = 1, . . . , Nl- Lo- 
cal effects also occur in short nanotubes, such as a sharp 
spike in the density of states (DOS), caused by defects 
at the caps. Such effects will be ignored here. Typically, 
the nanotube is assumed to be of infinite length, allowing 
continuous values of the wave vector along the nanotube 
axis. This causes possible wave vectors to lie in 'sub- 
bands'. The subbands can cut through the fermi points 
of graphene causing the tube to become metallic. This 
can be shown to be the case for nanotubes of the type 
(n, m), where n — m is a multiple of three. If this is not 
the case there is a nonzero band gap and the nanotube is 
semiconducting. Here we consider 'zigzag' semiconduct- 
ing nanotubes of the form (n, 0), with n not a multiple 
of three. The discrete wave vectors are then 

(18) 



fej = 



M 

with /i = —(n — 1), 0, 1, 2, n and 

7T 



-II 



< 



(19) 



The momentum vectors associated with these subbands 
are highlighted in Fig. [2] for a (7,0) nanotube. For zigzag 
nanotubes k± corresponds to k y and fc|| corresponds to 
k x . It should be noted that subbands /j, and — jj, both 
have the same energy and this degeneracy is referred to 
as the 'valley' degeneracy. Combined with the two elec- 
tron spins this gives a degeneracy of four for each energy 
value, except for /1 = and /i = n that only have spin 
degeneracy. The energy of the subbands is 



E± T (k ll ,»)=E ± (k ll 



K2 
IK, 



(20) 



These are plotted for a (7,0) nanotube in Fig. [31 which 
has a bandgap of 1.43eV. 




FIG. 3. (color online) The dispersion relation for a (7,0) nan- 
otube. The possible transitions caused by light linearly po- 
larized parallel to the nanotube are shown by arrows. Here 
the bandgap is 1.43eV, when 70 is taken as 2.89eV. 



III. THE OPTICAL ABSORPTION OF CARBON 
NANOTUBES 

The Hamiltonian of a nanotube interacting with an 
electromagnetic field is H = Hq + Hp + Hi, with 



H, 



e 



Ap, 



(21) 



being the interaction term and Hp representing the field 
Hamiltonian. Here, we define e as the magnitude of the 
electron charge and are using SI units. Each field mode is 
characterized by its angular frequency uj and further pa- 
rameters, which define the mode's polarization and prop- 
agation direction. The field vector potential operator is 



A = 



cM A+c 



+ — \UJt 



A-e iwt ) 



(22) 



We will consider the initial and final state of field to be 
a coherent monochromatic state |a Wo )> with an angular 
frequency of loq and a mean photon flux of F photons per 
unit time. This state satisfies the equation a w |a Wo ) = 
a\a Uo ), with a = v2ttFS(uj — ujq) and being the field 
mode's destruction operator—. This allows us to give 



(o^JAIa^) 



JUJ t 



(23) 
(24) 



Using time-dependent perturbation theory we find that, 
after time t, the initial state of the nanotube and field, 
I I oz-uo): is in the state |^I' /C ) \a Uo ) with probability 



!* ,c |(— A H 



P 



S(E' - E - hu Q ), (25) 



which is Fermi's Golden Rule. The transition rate for 
each electron in the state with energy E to each state 
with energy E' can be expressed as 



2tt 

IT 



—ihG 



6(E' - E -7kj ). 



(26) 
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To calculate the optical absorption of a carbon nan- 
otube of length L, the interaction term, ihG = 
iri(*' c |A+.V|* t '), needs to be found with spatially 
changing field. Here we are assuming that the state is 
coherent over the entire length of the nanotube. To cal- 
culate G we define 



v A (k)^c ik( >q, 
q 

v B (k) = -Eo-% 



(27) 
(28) 



with q summing over the three vectors pointing from an 
atom in the A sublattice to its neighboring three B lattice 
atoms. We will furthermore use the matrix element, M = 
(p z (r)|V z \p z (v — q z ), with q z being the vector between 
two neighboring atoms such that the z-axis is aligned 
along q 2 . The value we will later use for this is given 
by^ 



M = 2aj Q m e /h 2 V3. 



(29) 



Each of the unit cells in the nanotube extends over a dis- 
tance of |T| w 0.43 nm along the nanotube and approx- 
imately a nanometre across. This is much smaller than 
the light's wavelength and spatial variations. Therefore, 
we can assume that the electromagnetic field is constant 
across each of the nanotube's unit cells. There are Nl 
of these unit cells along the nanotube's length and in 
each one, labeled by an integer I. the field is given by 
A+ = A+(Z|T|-L/2). 

An expression for G can then be calculated and sim- 
plified into the form (see Appendix [A| 



G = — D(k',k). 



-D(k',k). 



E 

L/2 



^i( a ; v / 3-L/2)(fe||-fcf | ) A + 



;=-l/2 



d/e i'(fcii-Mi)A+(0 



(30) 
(31) 



where 



3=1 

(1+c -ia(k'-k).(y3/2,l/2) )v A (k) 

- c^(k')c^(k)e- ij ' a( ^- fei) e- m(fc !r fc n)/^ 

(1 + c --(k'-k).(y3 A l/2) )v A (k) * ); (32) 



and D x<y are given in Appendix [AJ 

This result coincides with that of RefJ^ when is 
the same for all I. In Eq. (f3"T|) the D gives the selection 
rules for possible transitions between bands // and /x. In 
particular D z (k' , k) is negligible if // ^ /j, and for the 
other components of D(k',k) to contribute we require 
p! = /i± 1. For a uniform field across the nanotube these 
give the possible transitions for a field polarized paral- 
lel and perpendicular to the nanotube, respectively. In 
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FIG. 4. (color online) A contour plot of D z (k, k) for graphene 
is shown together with the nanotube's subband lines. It is 
given in terms of the constant M from Eq. (1291 1, 
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FIG. 5. (color online) a) D x values for the different transitions 
that can occur when the electric field is perpendicular to the 
nanotube. M is given by Eq. (|29[) . b) D z for transitions 
allowed with an electric field parallel to the nanotube. 



Fig. 2] we have plotted D z (k, k), and in Fig. [5] D x ^ z is 
plotted. These expressions correspond to direct transi- 
tions, i.e with = fc||, which is an approximation of 
momentum conservation and will be discussed later in 
this section. 



Although the values of D x and D y show a transition, 
the induced local field creates a depolarization effect^— 
that reduces D x and D y to give a negligible contribution 
to the absorption. This allows us to focus on the D z term 
and simplify G to 



G«-j-£ 2 (k',k) 



A+.dre is(fc n 



(33) 



with s denoting the length along the nanotube. This also 
restricts the transitions to those with \i! = /.t. 

It is the line integral in Eq. (|33|) that is responsible 
for momentum conservation in the system. The photon 
momentum is much smaller than the crystal momentum 
and typically only direct transitions are considered, i.e. 
fc| | « Here we will make this assumption, however the 
change in momentum can not be completely neglected 
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since any change can make a major difference to the line 
integral in Eq. ([55]) . This is especially true when the field 
oscillates along the nanotube. The energy of a direct 
transition is given by E g (k) = E + (k) — E~(k). Since 
L> z (k',k) rj D z (k, k), when k',, « fey, we will make this 
substitution and further simplify D z (k, k) = D z (k) to 
give 



D z (k) 



1 0k 



(34) 



We define j4jj"(s)cis = A + .dr to be the field potential 
along the nanotube and use the discrete momentum val- 
ues, fc|| = 2irj/L and k'^ = 2irj / /L, with integers j and 
j'. The line integral can then be expressed as 



S(k\\ - fci. 



(1/L) ^A + .dre is(fc '^ fc !i) 



L/2 



(1/L) / dsA+(s)c is(fc "~ fc i|) 

' -L/2 



(35) 
(36) 



L/2 



(1/L) / dsA+(s)c i2 ™ ( ^" )/L . (37) 

J-L/2 



This expression is simply the coefficient in the Fourier 
series for Au~(s). Since the photon momentum is very 
small in comparison to the crystal momentum the only 
relevant coefficients will have very small values of j' — j 
relative to Nl- Every electron transition in the nan- 
otube then needs to be considered to calculate the over- 
all absorption rate. This leads to a length dependence 
on the absorption. We initially consider discrete states 
and corresponding ku values. The transition rate given 
by Eq. (|26p is summed over all possible initial and final 
states to give 



—EE EE^w-wi' 



S(fc||-fcf|) 8(E g Qt)-hua 



(38) 



In this equation di refer to the degeneracy of the initial 
state. For any value of k\\ the sum over causes k\ \ — fen 
to take all of the low values that are relevant. This sum 
is also independent of k\\ and allows us to define S = 

J2j \S(2nj / L)\ 2 , which can be rewritten using Parseval's 
theorem to be 



5= (1/L) 



,L/2 




/ ds 




J-L/2 



(39) 



The total absorption rate is then 

™ n t 2 

E E^H^wi 

di fj,=— n+1 fen e 



SS(E g (k) - hu ) 

[L/2 

« U(ljq) / ds 

J-L/2 



AUs) 



(40) 
(41) 



where we define 



he 2 



C/M = E E dk^\D z (k)\ 2 6(E g (k)-hu;o). 

^ p=-n+l e 

(42) 

The field was defined to be a coherent state with a pho- 
ton flux given by F photons per unit time. We divide the 
transition rate, given in Eq. (|4"Tj) . by this flux to obtain 
an estimate for the probability of one photon exciting a 
single electron. This expression gives a probability that 
increases linearly with nanotube length. This is certainly 
suitable up to the coherence length, L c , however not for 
long nanotubcs. So far we have considered the length of 
the nanotube to be smaller than the coherence length. 
For long nanotubes we can consider the whole nanotube 
to be composed of coherent segments. This leads to an 
exponential increase in the absorption with the nanotube 
length. Here, to calculate this quantity we find the prob- 
ability of not exciting any electrons, which is the product 
of (1 — u>l c /E) w e~ WL "l F for each segment. Hence, the 
probability of exciting a single electron, in a nanotube 
of length L, with each photon can be estimated by the 
expression 



( -UM f L/2 , 

7] = 1 — cxp | — — / as 



F 



L/2 



4(5) 



(43) 



Note that this expression is actually independent of the 
coherence length. 

Broadening effects can be included in this by substi- 
tuting the Dirac delta function, from Eq. (|4"2")l with a 
Lorentzian function, 



5{E g - hoj) 



7r((L 9 -M 2 +r 2 ) 



(44) 



that has a broadening parameter, T. This parameter can 
include the broadening due to multiple effects, includ- 
ing the electronic state's decay. In carbon nanotubcs 
the state decay occurs on a picosecond time scaled. If 
we take a range of O.lps to 2ps the required broadening 
ranges from T = O.OlcV to T = O.OOleV. In order to com- 
pare our results with previous wor k 25 ' 51 we will choose 
in the following to use a parameter of T = 0.0 leV. 

So far we have assumed the light to be in a pure state 
consisting of one specific wavelength. We expect to have 
a range of wavelengths present and to deal with this we 
assume the light is in a probabilistic mixture of coherent 
beams, each with a photon flux F. These are weighted 
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by a lineshape g(uj), satisfying J dujg(uj) = 1. The light's 
state is then J duig(u!)\a ul )(a ul \ and the expected absorp- 
tion probability is 



(45) 



r\ = / du)g(<jj)r). 



In the following we take g to be a uniform lineshape be- 
tween two energy values. This is equivalent to taking rj 
to be the average transition probability over a range of 
energies. 



IV. OPTICAL NANOFIBER PHOTON 
ABSORPTION INTO A CARBON NANOTUBE 

We now extend the calculation for the absorption 
around optical fibers and particularly nanofiber o 52 ' 53 . A 
review of these subwavelength diameter waveguides can 
be found in Rcfs^~— . They are made of a silica core and 
have diameters as small as 50 nm. For fibers of this size 
a high proportion of the light field exists outside of the 
fiber's core. This means the field is easily accessible and 
we consider positioning the carbon nanotube near the 
nanofiber. The use of fibers allows the interaction to be 
enhanced due to the transverse confinement of the field. 
Altering the nanofibers properties also allows us to tailor 
the field. Fibers with a smaller diameter have a larger 
evanescent field but also suffer from higher losses. We 
will consider a cylindrical nanofiber core with a radius of 
R and cladding provided by the vacuum, with refractive 
index n,2 = 1. The refractive index of the silica core is 
m = f .45 and the material absorption of the silica is neg- 
ligible over the short distances being considered. Silica 
core fibers with subwavelength diameter are single mode 
fibers, i.e. the only mode present is the HE fundamental 
mode (see Refj£& for a general single mode condition). In 
the following we will adopt a scheme used in Refs i 34 ' 46 
to quantize the field. The field potential operator for the 
nanofiber 's fundamental mode is then 



E 



hf3> 



(46) 



This is given in terms of cylindrical coordinates, with z 
being the coordinate along the fiber and ip the azimuthal 
angle. The light's angular frequency is cj. Its propaga- 
tion direction is labeled with / = ±1 and j3 refers to the 
longitudinal propagation constant. We find the value of (3 
by numerically solving the fiber eigenvalue equation (see 
Eq. (|BT]) in Appendix [Bj. The derivative in Eq. (|35]l. 
P' , is taken with respect to w, and a m are the photon 
annihilation operators, with m = {uj,f,p) characterizing 
the separate modes. Furthermore, e m are the electric 
field profiles of the guided mode that can be found by 
solving Maxwell's equation a 34 ' 57 and A gives a normal- 
ization factor. The expressions for the mode profiles and 
A are given in Appcndix[B] The polarization can be right 
or left circular labeled by p = ±1. For a single mode 
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FIG. 6. (Color online) The electric field, E, of the circularly 
polarized HE mode taken for a constant z, to give a cross- 
section of the fiber. The contour gives E z and the arrows 
represent the x and y components. The longer the arrow the 
stronger the field. The field has been divided by the con- 
stant Eq = yj Fuiohf} 1 1 (2e$A). In this particular case we took 
the fiber diameter to be 250nm and the light's wavelength as 
868nm. 



of monochromatic coherent light with m = (ujQ,f,p) we 
have 



l (r, <p)e 



i(fPz+ P <p) 



(47) 



Fig. |6] provides a slice of the classical field at one instant 
in time. This field can be seen to extend far away from 
the nanofiber and vary dramatically with position. 

This contrasts with the simpler case, that has previ- 
ously been studied, of a plane linearly polarized light 
beam that is given by 



.+ _ 



AiruieocA 



(48) 



across the whole nanotube, where the beam has a finite 
cross-sectional area of A' . This gives 



Fh 



2ujq6qcA' 



(49) 



For fibers larger than 100 nm in diameter the photon 
losses are small and can be ignored over short distances. 
In our calculations we will use nanofibers of diameter 
250 nm. Furthermore, we will focus on the forward prop- 
agation and right polarized guided modes, i.e. / = p = 
+1. All other modes are related to our results by sym- 
metry. The value of G is then highly dependent on the 
way the nanotube is orientated relative to the nanofiber 
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v[%] 

a) 0.14 



FIG. 7. (color online) Possible orientations of a straight nan- 
otube relative to a fiber of radius R. <j> labels the angle be- 
tween the nanofiber and nanotube. 



and can be calculated to be 



G 



Fh/3' 



2u) e A 

* c (k')|e m citv i(ff}z+pv) .V|*"(k)). 



(50) 



Since the field strength drops off exponentially the high- 
est value for the coupling will be achieved by having the 
nanotube as close as possible to the fiber. In our exam- 
ples, the distance between the nanotube's center and the 
surface of the nanofiber is chosen to be 1.25 nm. The nan- 
otubes we consider always have a radius less than 1 nm 
so this distance avoids any contact. There are two ori- 
entations we will consider. The first is that of a straight 
nanotube, of length L, oriented at an angle <f> relative to 
the nanofiber which includes parallel and perpendicular 
orientations as illustrated in Fig. [Jj 

For 2/xm nanotubes perpendicular to the fiber, the ab- 
sorption probabilities as defined by Eq. (|43[) for different 
wavelengths of light and different zigzag nanotubes are 
shown in Fig. HI We do not consider the absorption of 
photons with energies greater than 6eV since these are 
not visible and require the addition of the higher energy 
a orbitals for accurate results. Distinct absorption peaks 
are clearly visible and the largest absorption occurs for 
a (11,0) nanotube. The absorption for a nanotube in a 
linearly polarized coherent plane wave [see Eq. (|48|)] is 
also shown in Fig. [8] This beam has a cross-sectional 
area of 4fim 2 and exhibits the same absorption peaks 
as the fiber, but varies less with the light's frequency. 
It can be seen that the (11,0) nanotube has its small- 
est energy transition dramatically reduced. This extra 
effect is caused due to larger evanescent fields, for an 
increasing wavelength relative to the fiber radius. This 
reduces the field intensity and absorption. The quan- 
tum efficiencies are a similar order of magnitude as those 
observed experimentally for plane wave o 16 ' 49 ' 58 . The dif- 
ferent nanotubes show shifted absorption peaks. These 
can be further adjusted with external fields or choosing 




400 600 800 1000 1200 1400 
Wavelength [nm] 

FIG. 8. Photon absorption probabilities, 77, for a 2 /im long 
nanotube, perpendicular to the fiber, at different photon en- 
ergies. The nanotubes considered are a) (7,0) b) (8,0) and 
c) (11,0). In each case the solid lines refer to the absorp- 
tion of circularly polarized light guided by the nanofiber and 
the dashed lines represent the absorption for a plane coherent 
light beam (without a fiber) that is linearly polarized along 
the nanotube. The smallest and second smallest transitions, 
En and E22, are indicated. We have used a broadening pa- 
rameter of T = 0.01 eV. 



other nanotubes^—. The resonant energy values are un- 
changed with the orientation and this allows us to choose 
a range to average over as a general measure of absorp- 
tion. We chose to calculate the mean absorption fj over 
the (7,0) nanotube's lowest absorption energy, particu- 
larly we chose a range of leV from 1.3cV (953nm) to 
2.3eV (539nm). The resulting fj is approximately inde- 
pendent of T in a range of T = O.OlcV to T = O.OOleV 
deviating only by a few percent. 

The corresponding mean absorption against nanotube 
length, for the lowest energy transition, is shown in Fig. [5] 
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FIG. 9. Average photon absorption probability, rj, for a 
straight (7,0) nanotube of length L. At L = from top 

to bottom, the angles between nanotube and fiber are <f> = 
(parallel), —7r/32, —37r/8, 7r/2 (perpendicular), tv/32, n/8. 
The mean absorption has taken over a leV region, from 1.3eV 
to 2.3eV. 



for various angles between the straight nanotube and 
nanofibcr. The results show that the absorption con- 
verges to a maximum value as the length is increased, 
unless the nanotube is parallel to the fiber. The nanotube 
perpendicular to the fiber has a very strong absorption 
for short lengths. In this situation we see the absorption 
increasing strongly with nanotube length which is due to 
the linear increase in electron number. As the length in- 
creases further this effect is counterbalanced by the fact 
that the field strength decreases exponentially away from 
the nanofiber. The absorption hardly increases at all af- 
ter the nanotube exceeds approximately 2 fj,m. However, 
over these short distances the absorption of the perpen- 
dicular nanotube can be improved upon by shifting the 
nanotube slightly away from a perfectly perpendicular 
arrangement. The parallel orientation increases slowly 
but docs not peak. This effect will be discussed later 
in this section and wc will find that the probability can 
be enhanced by spiralling the nanotube to combine both 
effects. If linear polarized light was used instead of circu- 
lar polarized light the absorption could be twice as high 
depending on the nanotube's position in the nanofibcr 
plane. We also see a difference between angles of ±7r/32, 
with the higher absorption being dependent on the light 's 
polarization and propagation. 

The strong absorption for a perpendicular nanotube 
is limited by the drop-off in field strength. However, 
this can be prevented by maintaining a constant dis- 
tance between the nanotube and nanofiber center, R n . 
The nanotube can locally approximate a perpendicular 
nanotube by spiralling around the nanofiber, as illus- 
trated in Fig. [TUl Although this bending does alter the 
electronic and optical properties these effects are small 
and can be safely ignored tiered. We define a 'wind- 
ing number', W , for the spiral as the number of loops 
per unit length along the z axis. This winding num- 
ber is equal to W = 1/di where di is the z-distance for 
one loop. An angle is also formed between the spiralling 




FIG. 10. (Color online) Nanotube spiralling around a 
nanofiber. The length of one loop, along the z axis, is la- 
beled as d(. 




FIG. 11. Average photon absorption probability, rj, for a (7,0) 
nanotube of length L coiled around the fiber. The average is 
taken from 1.3eV to 2.3eV. From top to bottom the winding 
numbers are -0.0016nm _1 , -0.0008 nm~\ nnT 1 (Parallel), 
0.0016 nm _1 and 0.0008 nm" 1 . 



nanotube and the nanofiber's direction, which is given by 
<t> s = arctan(27rVFi?„). Since these spiralling nanotubes 
can interact with the field over an arbitrary length their 
absorption's approach 100% given sufficient length and 
an allowed transition. 

The average absorption probabilities in this case are 
plotted in Fig. [TT] and show a steady increase in the ab- 
sorption probability with nanotube length. The paral- 
lel nanotube is also shown. This demonstrates that the 
spiralling nanotubes can have higher absorption proba- 
bilities than the parallel configuration. In Fig. [T^] we 
have plotted the average absorption probability against 
$ s for nanotubes of different lengths. An optimal spi- 
ralling rate to enhance the absorption can be seen. Wc 
found that the optimal value of this winding rate is 
W op t = e™ /(2irRne™). This was obtained by maximizing 
the alignment between the nanotube and e m . 
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FIG. 12. Average photon absorption probability, rj, for (7,0) 
nanotubes coiled around the fiber for different winding pa- 
rameters. From top to bottom the nanotubes have lengths 
lOOOOnm, 5000nm and lOOOnm. The average is taken from 
1.3eV to 2.3eV. 



APPLICATIONS 




FIG. 13. (Color online) Illustration of a possible photodc- 
tector. Here a 'forest' of aligned semiconducting nanotubes 
(thick dark lines) are grown between two electrodes. The 
nanofiber is positioned between these electrodes. Once a light 
field excites an electron the resistance between the electrodes 
drops dramatically, which allows the photon to be recorded. 



The nanotube-nanofiber setups discussed in the pre- 
vious sections open up possibilities for a range of ap- 
plications, particularly highly sensitive photodetectors. 
These systems would detect light guided within a fiber. 
In this section we discuss the possibilities. Note that we 
have only considered the quantum efficiency of the ab- 
sorption and that detection of the charge excitations is 
beyond the scope of this paper. However, certain nan- 
otube properties such as ballistic electron transport and 
low capacitance should be a great advantage for this de- 
tection. The bandgap of carbon nanotubes decreases 
with the nanotube size, so for optical and near infrared 
wavelengths small diameter nanotubes are required. This 
rules out the possibility of encasing a nanofiber within 
a nanotube. Instead, a practical setup is given by ar- 
ranging N horizontal nanotubes in a parallel array and 
placing the nanofiber orthogonally on top of the array. 
Based on current nanotube arrays, the density of nan- 
otubes would be 1 — 100 nanotubes per fim £2r£l. We 
will use rj as a measure of the absorption probability for 
one nanotube. The photon absorption probability of each 
nanotube is then, in case of a (7,0) nanotube, given by 
the top line in Fig. [9] and the overall absorption proba- 
bility is 



Vtot = 1 - (1 - v) 



N 



(51) 



Taking rj = 0.00015 (see Fig.© this leads to r] to t > 95% 
for N > 20000, a value greatly exceeding those of cur- 
rently available APDs^. For N > 40000 the efficiency 
exceeds 99% which can currently only be achieved by 
highly complex superconducting detectors^. The ad- 
vantage of our setup is that it can be operated at room 
temperature. Each nanotube would require a length of 
2/xm and has to be connected at the ends by electrodes^ 
which collect the excited electrons via an applied volt- 



age. Although this should be possible in the near future, 
current technology cannot generate an array of unique 
nanotubes. 

Aligned vertical nanotubes can also be grown on a con- 
ducting substrate, that can then serve as one electrode. 
The nanofiber can then be placed orthogonally to the 
nanotubes and the remaining ends of the nanotubes con- 
nected to an additional electrode (see Fig. [T3|). The di- 
ameter of the nanotubes in this case can be in the range 
of 1 ± 0.5 n m 67 i 68 . Recently, progress has been made in 
the generation of such semiconducting nanotube 'forests', 
although a semiconducting nanotube purity of 100% has 
yet to be achieved reliably£2r2I. These nanotube systems 
typically have a density of 10—10000 nanotubes per /mi 2 
— ^-i. The nanotubes are distributed uniformly over a 
selected region and we assume that they are a uniform 
mix of semiconducting zigzag nanotubes with a diame- 
ter in the range 1 ± 0.5 nm. We calculated the overall 
absorption probability when we have a forest that ex- 
tends a distance of 500nm from the nanofiber and 15 /im 
along its length, with a density of 900 nanotubes per 
/im 2 . The results with nanofibcrs that have diameters of 
250 nm and 400 nm are shown in Fig. Q31 These absorb 
light of a wide range of wavelengths, that can be selected 
by the nanotubes present and choice of nanofiber diam- 
eter. A typical absorption probability of r) to t > 50% can 
be seen, for 250 nm diameter fibers, and by extending the 
system's length from 15 /im to 50 /im this is increased to 
rjtot > 95%. Nanotubes around a 400 nm fiber are also 
seen to absorb light at wavelengths that are typically 
used for optical communication. Due to the nanotube's 
bandgap dependence on external fields there is also the 
possibility of adjusting the absorption frequencies by us- 
ing an external field. 

An additional possible setup is given by arranging the 
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FIG. 14. Absorption probability for a circularly polarized 
photon in a nanofiber laid inside a 'forest' of nanotubes (see 
Fig. I13[) . This is given for 2 /im long vertically aligned nan- 
otubes. The nanotubes are in a region that extends 500nm 
away from the fiber and 15/tm along its length. The density 
of the array is taken as 900 nanotubes per /im 2 . The fibers 
diameter is taken to be either 250 nm (black line) and 400 nm 
(gray line). A broadening parameter of F = 0.01 eV was used. 



have found a simple expression for the absorption prob- 
ability, which is independent of the coherence length. 
Currently, the coherence lengths of carbon nanotubes is 
an area of extensive research with results ranging from 
10 nm, at room temperature, to a few microns^—. They 
seem to be highly dependent on the temperature, im- 
purities, defects and surrounding fields. Once excited, 
the radiative lifetimes of the excitations have been ob- 
served to range from 3 to 100 n o 28 i 76 . The nonradiative 
decay seems to be much faster, of the order of a few 
picoseconds^. Introducing excitons into the model and 
analyzing the dynamics of quantized single photon states 
within a carbon nanotube is an interesting area of further 
research. Yet, the results here are an important step to- 
wards calculating the absorptions within nanostructures 
and are of importance to future nanotube optoelectonics. 

The research of UD was supported by the Centre 
for Quantum Technologies, National University of Sin- 
gapore and the ESF via EuroQUAM (EPSRC Grant 
No. EP/E041612/1). SB acknowledges support from the 
EPSRC Doctoral Training Accounts. 



nanofiber and the nanotube parallel to each other. Tak- 
ing 100 nanotubes of length L = 1 mm parallel to the 
fiber and using rj = 0.07 (sec Fig. IT!]) we obtain an overall 
absorption probability of r\tot > 99% which again greatly 
exceeds that of standard APDs. 

As a final setup we consider the coil geometry shown 
in Fig. [10] which has a high absorption probability of 
up to 100% for long nanotubes. However, producing 
such a setup in a laboratory is rather challenging with 
current technology. This setup also allows for further 
specification of the absorbed light's polarization or prop- 
agation direction with the choice of winding number. 
The winding also dramatically reduces the length of the 
system. For a nanotube with a winding number of 
W = — 0.1 nm -1 , the average absorption between 1.5eV 
and 2.5eV exceeds 50% when the nanotube's length is 
5 mm. For the 250 nm diameter fiber this only extends 
64 fan along the fiber. 



VI. SUMMARY 

In this paper we have calculated the probability of ab- 
sorbing a photon with zigzag carbon nanotubes. The 
light field is allowed to vary along the nanotube, i.e. no 
dipolc approximation is made, which has enabled us to 
treat the absorption of light from optical nanofibcrs. We 
found that there is a strong dependence on the system's 
geometry and have devised setups for high absorption. If 
we spiral the nanotube around the fiber, we find that an 
absorption of circularly polarized light, arbitrarily close 
to 100% can be achieved. For straight nanotubes, that 
are not parallel to the fiber, we find the absorption prob- 
ability converges as the nanotube's length increases. We 



Appendix A: Calculating the Optical Matrix 
Element 

The matrix element for the interaction between can be 
found by substituting in the wavefunctions to give 

G = (* c (k')|A + .V|**(k)} 
= ]T cr(k / K(k)(pJ(k / )|A+.V|p t ,(k)) 



s,t=A,B 

] 

Ncells 



±- £ c r(k'K(k) 

:,t=A,B 

J- e ikr2 - ik ' ri (p z (r - n)|A+.VMr - r 2 )) 

ri6R 3 ,r 2 GR t 
1 

= N G N L 
cl(k'K(k)e 



E 

riERA,r 2 €R B 



(p z (r-n)|A + .V|^(r-r 2 )) 
c£(k')c^(k)e ik - ri - ik '- r2 

(p z (r-r 2 )|A+.VMr-rx)) 

' aN G N L 



(cX(k')c^(k) 



e ^ i(k '- k) - ri A+(r 1 ).v j4 (k) 



i-iGRa 

+ c£(k')c^(k) 

c - i ( k '- k )- r2 A+(r 2 ).v B (k)). (Al) 

r 2 <£B B 

Here we have assumed that the orbitals are symmetric 
and that A + is constant across each of the nanotube's 
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unit cells. The expression for G can then be split into 
separate unit cells and directions 



G — G x + G y + G z 

G z = Mv 5 (^ c i(alV *- (L/2))(k v- k 'n ) 



(A2) 



2anN T 



1=1 



At(laV3-(L/2))) 

n 

^(cX(k')4(k)e- ija ^-^) 

(1 + e -ia(k'-k).(V3/2,l/2) )t; A (k) 

-c% (k')c^(k)e- i3Q( ^- fci) e ia(fc ii~ fe !i)/v / 3 

(1 + c _ ia(k '_ k) . ( y3Al/2) )v A (kr) 
1 w i 



1=1 



A+(W3-(L/2))). 



(A3) 



G rf = ^ c i(fcn-fc(|)»aA/3-(L/2)) 



2anN L 



1=1 



A+(laV3-(L/2))) 

n 

E c A(k')^(k)« 27rj y n) (k)e- ya ( fc -- fe ^+ 



d(27r0'+l/2)/«) 

c£(k')cl(k)(« 



(k)c _ija(fc ^ _fei) e~ iak '- (1/(2Vl)a/2) ) 



d(2nj/n) 



(k)e 



— ij'o^^L — iofej|/ v3 



^(2 ff0 - + i/2)/ n )(k)e-^^-^)e-- k '( 5 /( 2 V3)-V 2 )) 
— D d (^e i(fc|l " fc !i )(Zav/I_(L/2)) 

A+(W5-(L/2))). 



i=i 



(A10) 



Hence, we obtain 



G = ^ c i(fc N-Mi)(^-( i / 2 ))A+(/aV3 - (i/2)). 

(All) 



Z=l 



In order to calculate G x and G y we must take the cur- 
vature of the nanotube into consideration. To do this we 
use the method from Rcf^ and introduce the parameters 



Appendix B: Classical Nanofiber Field Modes 

For light of wavelength, A, and k = 2tt/X, the field 
parameters must satisfy the fiber eigenvalue equation^ 



>k.(-l/(2\/3), -1/2) ^270/271 



1) 



_j_ c iak.-l/(2 v / 3,l/2)(' c ±27ri/2)i _ jN 
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„,,l/2)(- e ±27ri/2,i _ ^ 



^ (e) (k) = i?, 
< W (k) = 



< W (k) = i^ 
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+ e" 


-ie v A 
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e [9 v^° 
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■i6 v B 
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e ie v*° 


+ e~ 






2i 




e i9 v+° 
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2i 



From these we calculate Gd= x y to be 



(A4) 

(A5) 
(A6) 

(A7) 

(A8) 

(A9) 



MhR) 
hRJiQiR) 



K[(qR) 



1 



2n\ qRKi(qR) h 2 R 2 



n{ - n\ K 



'MR) V 



2n 2 qRKi(qR) J 



P 2 ( 1 



n 2 fc 2 V<7 2 ^ 2 ^ 2 # 2 



]Va, (Bl) 



with J v referring to the Besscl function of the first kind 
and K v being the modified Besscl function of the second 
kind. By numerically solving Eq. (|B1[) the value of the 
propagation constant, j3, is determined. We also define 
the parameters h = (n\k 2 - /3 2 ) 1/2 , q = (/3 2 - n^k 2 ) 1 / 2 
and 



1 1 \ / ( J'l(hR) K[(qR) \ 

" \q 2 R 2 h 2 R 2 J ' \hRJi(hR) qRK\(qR) ) ' 

(B2) 

Inside the fiber (0 < r < R) the guided mode, m = 



(f,p), have the form and outside r > R 
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e™ = i[(l - g)K (qr) + (1 + g)K 2 (qr)}, (B6) 
= -p[(l - ff)Jifo(gr) - (1 + fl)Jf 3 (gr)], (B7) 



e " 1 = i |f^ [(1 -- 9)Jo(/ir) - (1 + ' 9)J2(H] ' (B3) e™ = &K 1 ( q r), (B8) 

( ) 

e ^ = ~~ P HJi(hR) ~ &'^°V ir ' ^ ~^ fl0^2(7w)]) These are normalized with a factor given by 

(B4) 

/>27T /«.R /»27T pOO 

m 2qKi(qR) / / n\\e\ 2 rdrdp + / n| |e| 2 rdrdtp = A. 

e z = f n t fh m Jl ( hr )> ( B5 ) Jo J ° Jo Jr 

pJi(hR) (B9) 
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